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1. Introduction 

The evaluation of the cross-sections for hadron-hadron iteractions needs the sufficiently pre¬ 
cise knowledge of parton distribution functions (PDFs) and parton fragmentation functions (FFs), 
which are the important part of any cross-section. The properties of PDFs and FFs can be taken 
from processes of the deep-inelatic scattering (DIS) and e + e~ -collisions, respectively. In this re¬ 
port we will concentrate only for the high-energy limits of PDFs and FFs, which are needed for 
modern experiments studied on LHC collider. 


1.1 PDFs 

The experimental data from HERA on the DIS structure function (SF) Fj [0H0]- its deriva¬ 
tive dlnF2/<9ln(l/x) [4] and the heavy quark parts F™ and F^ h [||] enable us to enter into a very 
interesting kinematical range for testing the theoretical ideas on the behavior of quarks and gluons 
carrying a very low fraction of momentum of the proton, the so-called small-* region. In this limit 
one expects that the conventional treatment based on the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equations [j6|] does not account for contributions to the cross section which are 
leading in a s ln(I/*) and, moreover, the parton densities are becoming large and need to develop a 
high density formulation of QCD. However, the reasonable agreement between HERA data and the 
next-to-leading-order (NLO) and next-to-next-to-leading-order (NNLO) approximations of pertur¬ 
bative QCD has been observed for Q 2 >2 GeV 2 (see reviews in [7j| and references therein) and, 
thus, perturbative QCD could describe the evolution of F 2 and its derivatives up to very low Q 2 
values, traditionally explained by soft processes. 

The standard program to study the * behavior of quarks and gluons is earned out by compar¬ 
ison of data with the numerical solution of the DGLAP equation [B] 1 by fitting the parameters of 


the PDF A-profilc at some initial <2o and the QCD energy scale A [fql-[|l 2[]. However, for analyzing 
exclusively the low-* region, there is the alternative of doing a simpler analysis by using some of 


the existing analytical solutions of DGLAP evolution in the low-* limit [ 13]—[ 16]. This was done 


so in [ 13] where it was pointed out that the HERA small-* data can be interpreted in terms of the 


so-called doubled asymptotic scaling (DAS) phenomenon related to the asymptotic behavior of the 
DGLAP evolution discovered many years ago Jl7t ]. 

The study of [ 13 ] was extended in [ |T-i| . [ffij [i~6| | to include the finite parts of anomalous dimen¬ 
sions of Wilson operators 2 . This has led to predictions [ L5, 1_6] of the small-* asymptotic PDF 


form in the framework of the DGLAP dynamics starting at some Qfi with the flat function 


fa(Qo) = A a (hereafter a = q,g), 


( 1 . 1 ) 


where f a are the parton distributions multiplied by * and A a are unknown parameters to be deter¬ 
mined from the data. 


'At small * there is another approach based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [|i|], whose 
application is out of the scope of this work. 

2 In the standard DAS approximation 0 only the singular parts of the anomalous dimensions were used. 
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1.2 FFs and average jet multiplicities 


Collisions of particles and nuclei at high energies usually produce many hadrons and their 
production is a typical process where nonperturbative phenomena are involved. However, for par¬ 
ticular observables, this problem can be avoided. In particular, the counting of hadrons in a jet 
that is initiated at a certain scale Q belongs to this class of observables. Hence, if the scale Q is 
large enough, this would in principle allow perturbative QCD to be predictive without the need to 
consider phenomenological models of hadronization. Nevertheless, such processes are dominated 
by soft-gluon emissions, and it is a well-known fact that, in such kinematic regions of phase space, 
fixed-order perturbation theory fails, rendering the usage of resummation techniques indispens¬ 
able. As we shall see, the computation of average jet multiplicities (AJMs) indeed requires small-* 
resummation, as was already realized a long time ago [|T^] . In Ref. |[T8||. it was shown that the 
singularities for * ~ 0, which are encoded in large logarithms of the kind I fx\n k { I /*), spoil per¬ 
turbation theory, and also render integral observables in * ill-defined, disappear after resummation. 
Usually, resummation includes the singularities from all orders according to a certain logarithmic 
accuracy, for which it restores perturbation theory. 

Small-* resummation has recently been earned out for timelike splitting fuctions in the MS 
factorization scheme, which is generally preferable to other schemes, yielding fully analytic ex¬ 
pressions. In a first step, the next-to-leading-logarithmic (NLL) level of accuracy has been reached 
[19, 20| [. In a second step, this has been pushed to the next-to-next-to-leading-logarithmic (NNLL), 
and partially even to the next-to-next-to-next-to-leading-logarithmic (N 3 LL), level [(2~jj| . Thanks to 
these results, we were able in [ [221 [u| l to analytically compute the NNLL contributions to the evolu¬ 
tions of the gluon and quark AJMs with normalization factors evaluated to NLO and approximately 
to next-to-next-to-next-to-order (N 3 LO) in the y/~a s expansion. The previous literature contains a 
NLL result on the small-* resummation of timelike splitting fuctions obtained in a massive-gluon 
scheme. Unfortunately, this is unsuitable for the combination with available fixed-order correc¬ 
tions, which are routinely evaluated in the MS scheme. A general discussion of the scheme choice 
and dependence in this context may be found in Refs. [|24|]. 

The gluon and quark AJMs, which we denote as (ni 1 (Q 2 )) g and ( nj r (Q 2 )) q , respectively, rep¬ 
resent the average numbers of hadrons in a jet initiated by a gluon or a quark at scale Q. In the 
past, analytic predictions were obtained by solving the equations for the generating functionals in 
the modified leading-logarithmic approximation (MLLA) in Ref. [25] through N 3 LO in the ex¬ 
pansion parameter i.e. through However, the theoretical prediction for the ratio 

r(Q 2 ) = ( n h(Q 2 ))g /( nh(Q 2 ))q given in Ref. J25| ] is about 10% higher than the experimental data at 
the scale of the Z° boson. An alternative approach was proposed in Ref. [|2fi||, where a differential 


equation for the gluon-to-quark AJM ratio was obtained in the MLLA within the framework of the 
colour-dipole model, and the constant of integration, which is supposed to encode nonperturbative 
contributions, was fitted to experimental data. A constant offset to the gluon and quark AJMs jet 


multiplicities was also introduced in Ref. [27]. 


Recently, we proposed a new formalism [28, 22, El] that solves the problem of the apparent 


good convergence of the perturbative series and does not require any ad-hoc offset, once the effects 
due to the mixing between quarks and gluons are fully included. Our result is a generalization 


of the result obtained in Ref. [25]. In our new approach, the nonperturbative informations to the 
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gluon-to-quark AJM ratio are encoded in the initial conditions of the evolution equations. 

This contribution is organized as follows. Section 2 contains general formulae for the Q 2 - 
evolution of PDFs and FFs. The generalized DAS approach is presented in Section 3. Sections 
4 and 5 contain basic formulae of <2 2 -dependence of FFs at low x and the AJMs, respectively. In 
Section 6 we compare our formulae with the experimental data for the DIS SF F 2 and the AJMs 
and present the obtained results. Some discussions can be found in the conclusions. The procedure 
of the diagonalization and its results for PDF and SF Mellin moments can be found in AppendixA. 


2. Approach 

Here we breifly touch on some points concerning theoretical part of our analysis. 

2.1 Strong coupling constant 

The strong coupling constant is determined from the renormalization group equation. More¬ 
over, the perturbative coupling constant a s {Q 2 ) is different at the leading-order (LO), NLO and 
NNLO approximations. Indeed, from the renormalization group equation we can obtain the fol¬ 
lowing equations for the coupling constant 


1 


1 


^°(2 2 ) a\°{M 2 

at the LO approximationm and 
1 1 


= w»[g 


^ LO (£ 2 ) a" LO (M 2 ) 


+ b \In 


^ 0 (e 2 )(i+M? L0 (Mf)) i / e 2 

flN LO (M|)( 1+ ftiflN LO ((2 2 ))J P ° \ M 2 


( 2 . 1 ) 


( 2 . 2 ) 


at the NLO approximation. 

At NNLO level af NLO (<2 2 ) = a .\(Q 2 ) is more complicated and it is given by 


1 


1 


a s (Q 2 ) a s (M z ) 


+ b\ In 


a s (Q 2 ) h+b\a s (M^) + b 2 a 2 {M^) 
a / M z)\ l+b l a s (Q 2 ) + b 2 a 2 (Q 2 ) 


b 2 

2 


+ (*2-2) • / = A)ln ( ^2 


(2.3) 


The expression for I looks: 


/ = 


2 ( b 1 +2b 2 a s (Q 2 ) bi+2b 2 a s {M%,) 

—= arctan-—-arctan-—-— 

s/A V s/A s/A 


1 


: In 


b\ +2b 2 a s (Q 2 ) ~V~A b\ +2b 2 a s (M z ) + s/^A 


, s/—A \_b\_+2b 2 a s (Q 2 ) + s/^A b\ -\-2b 2 a s (M^) — s /—A 
where A = Ab 2 — b 2 and b\ = are read off from the QCD /3-function: 


for / = 3,4,5;A > 0, 
for / = 6; A < 0, 


j8 (a s ) = - foci 2 - j8i a] - f3 2 a 4 s +... = -j3 0 a 2 1 +b\a s + b 2 a 2 + ... 


(2.4) 


where 


_ 2 38 „ 2857 5033 , 325 , 

*=. 1-3 f' ' Si = 102 - t / ’ lh = — —w f+ i 4 f -- 


with / being the number of active quark flavours. 


(2.5) 
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2.2 PDFs and DIS SF F 2 

The DIS SF can be represented as a sum of two terms: 

F 2 (x, Q 2 ) = F 2 ns (x, Q 2 ) + F 2 (x, Q 2 ), 


( 2 . 6 ) 


the nonsinglet (NS) and singlet (S) parts. At this point let’s introduce PDFs, the gluon distribu¬ 
tion function f g (x,Q 2 ) and the singlet and nonsinglet quark distribution functions fs(x,Q 2 ) and 
fNs(x,Q 2 ) 3 : 

/ 

fs(*,<2 2 ) = Y, { i( x ’Q 2 ) = v (fQ 2 )+s(x 1 q 2 ), 

i 

f NS (x,Q 2 ) = £c 2 ff (x,Q 2 ), t? s (x,Q 2 )=fi(x,Q 2 )-jf s (x,Q 2 ), (2-7) 

i J 

where V(x,Q 2 ) = u v (jc,<2 2 ) + d v (x,Q 2 ) is the distribution of valence quarks and S(x, Q 2 ) is a sum 
of sea parton distributions set equal to each other. 

There is a direct relation between SF moments M„(Q 2 ) and those of PDFs 

M j (n,Q 2 ) = f l dxx f, - 2 F{(x,Q 2 ), fj(n,Q 2 ) = j dxx n ~ 2 fj(x,Q 2 ), ( j=NS,S,G ), (2.8) 

JO JO 

which has the following form 


(Q 2 ) = C NS (n,d s (Q z )) -fNs(n,Q z ), 


K(Q 2 ) = e{Cs{nMQ 1 )).f s {n,Q 1 )+C g {n,d s (Q 1 ))-f G ^Q 2 ) ) , 


(2.9) 

( 2 . 10 ) 


with Cj(n,a s (Q “)) (j = NS,S,G ) are the Wilson coefficient functions. The constant e depends on 
weak and electromagnetic charges and is fixed for electromagnetic charges to 


1 f 

<= W 

J a 


( 2 . 11 ) 


Note that the NS and valence quark parts are negledgible at low x and, thus, F 2 (x,Q 2 ) = 
Fi(x,Q 2 ) and S{x,Q 2 ) = f s (n,Q 2 ) =i q {n,Q 2 ). 

2.3 Q 2 -dependence of SF moments 

The coefficient functions C q {n 1 a s {Q 2 )) = Cs(n,a s (Q 2 )) andC ? («.a v (0 2 )) are further expressed 
through the functions B q \n) and B g \n), respectively, which are known exactly [29, Q 4 

C a (n,a s {Q 2 )) = 1 + a s ■ B { C P(n) + a 2 ■ B [2 \n) + @(a 2 ), (a = ( q,g )) (2.12) 

where is the Kronckei' symbol. 


1 Un I i ke the standard case, here PDFs are multiplied by jc. 

4 For the integral and even complex n values, the coefficients B',P (n) and Z^lin) (a,b = q,g) can be obtained using 
the analytic continuation & 
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The (^-evolution of the PDF moments can be calculated within a framework of perturbative 
QCD (see e.g. [30, 32]). After diagonalization (see Appendix A), we see that the quark and gluon 
densities contain the so called ” + and ” — ’’-components 


fa{n,Q 2 ) =tf(n,Q 2 )+f a (n,Q 2 ) (,a = q,g ), 


(2.13) 


which in-turn evaluated already independently: 


tf(n,Q 2 ) =f±(n,Qo)- 


~ ffi(G 2 ) ' 

.a*(Go). 



•Rf(ffiG 2 ), 


where 



(2.14) 


(2.15) 


is the anomalous dimensions of the ” + and ” — ’’-components, which are obtained from the 
elements of the martix of the LO anomalous quark and gluon anomalous dimensions. 

At LO, the normalization coefficients f f(n, Qq) have the form 

(2.16) 


where 5 


*7(ffiGo) = fs(n,Ql) ■ a n + f g {n,Ql) ■ p ni fRi,Go) = Go) ■ (1 - <4)+ f <?(ffiGo) • £«, 
%(n,Ql) = i a {n,Ql)-r a {n,Ql) (2.17) 


_ Yqq\n)-7+\n) o _ Y { qg(n) _ Ym{n) 

y { - W - yf ] {n) ’ ri 0) („) - r f („) ’ '' (») - rf (»)' 

Above LO, the normalization factors f f(n, <2q) become to be 

fRn,G5) = £(n,Qg)- (l-fl^G^zi'iR-fl^G^zEw) 
+^(n,G6)-«,(Go) + 

where (see Appendix A, where Z± ± = —V±]_, Z^ ± a = —V^_ a , (i = 1,2)) 


(2.18) 


(2.19) 



5 To conrary [po| we replace a n by j8». Another expressions for the projectors a n , /!„ and e„ can be found in |[7~3||. 
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and 


(0) / \ (0) / \ 

= 4± q (n) = z¥ ± (n), Z% >g (n) = Z^ q (n) • ^ W ~ \ ” , (i=l,2) (2.21) 

7VW-7± («) 

Here y±±(n) and j±L{n) are the elemens of matrices of anomalous dimensions, which have been 
obtained after diagonalization from («) (a,b = q,g ) (latter taken in the exact form from [f4|]): 

7-- («) = Yqq (n) ■«« + 7« W ■ £« + Y$ (n) • A, + y® (n) • (1 - a n ), 

7-|(«) = 7-1 (n) ~ Yqg(n)^j , 7+j(«) = 7?? (n) + - 7—0), 


7+10) = tS-O)- ( 7 ^ 0 )-^7^0)) • 

The function H+ (n. Q 2 ) up to NNLO may be represented as 


( 2 . 22 ) 


Hain-Q 2 ) = l+a s (Q 2 )(z£l(n)-z£l a (n)j + a 2 (Q 2 ) (z^(n) -zg^n) j + 0 (a]{Q 2 )) , 

(2.23) 

where (see [Btl] and Appendix A) 


Z ±±(p) = Z ±±0) + E Z ±i,q Z i±,<P = Z ±U n )+L Z ±lq Z *,q’ 

i=± i=± 

Yqq ] 0) ~ 7+ 0) Q) 


r (.) (n) _ /± o, w - 

2.4 Fragmentation functions and their evolution 


( 2 ) 


(2.24) 


When one considers AJM observables, the basic equation is the one governing the evolution 
of FFs D a (x.p 2 ) for the gluon-quark-singlet system a = g,q. In Mellin space, it reads: 


D q (co,ii 2 ) \ _ ( P qq (a),a s ) P gq (co,a s ) \ ( D q (co,p 2 ) 


d^i 2 l D g ((o,fJ, 2 


Pqg(o).a s ) P gg (co,a s ) I \ Dg(co,fi ) I 


(2.25) 


where P a b(co,a s ), with a, b = g,q, are the timelike splitting functions, 6 (0 = n— 1, with n be¬ 
ing the standard Mellin moments with respect to x. The standard definition of the hadron AJMs 
corresponds to the first Mellin moment, with co = 0 (see, e.g., Ref. [35]): 


{n h {Q-)) a = 


dxx m D a (x, Q 2 ) 


= D a ((O = 0,Q 2 ), (a = g,q). 


(2.26) 


. to=0 


The timelike splitting functions P a b((D,a s ) in Eq. ([2.25) may be computed perturbatively in a s , 


P a ^(D,a s ) = f j a k s +l pM{(D). 
k =0 


(2.27) 


The functions P'^ (to) for k = 0. 1.2 in the MS scheme may be found in Refs. | [36| , |37| , |38| | through 
NNLO and in Refs. [ ]I9| , ^(i| , 21 ] with small-x resummation through NNLL accuracy. 

6 P a b = —Ya bl^' where yj. are the timelike anomalous dimensions. 
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2.5 Diagonalization of FFs 

As it was in the spacelike case (see subsection 2.3 and Appendix A), it is not in general possible 
to diagonalize Eq. ( 2.25 ) because the contributions to the timelike-splitting-function matrix do not 
commute at different orders. The usual approach is then to write a series expansion about the 
LO solution, which can in turn be diagonalized. One thus starts by choosing a basis in which the 
timelike-splitting-function matrix is diagonal at LO (see, e.g., Ref. [R(l| and Appendix A), 


P(co,a s ) = 


P ++ ((D,a s ) P~ + (a>,a s ) 
P + _(CD,a s 


*s ) P — ( CO , Cl s ) 


= a c 


'ifi(ffl) 0 


0 


P { °l(co) 


+ a 2 P W (co) + 0(a 3 s ), (2.28) 


with eigenvalues P'l+icD). In one important simplification of QCD, namely JP = 4 super Yang- 
Mills theory, this basis is actually more natural than the (g, q) basis because the diagonal splitting 
functions P±± (ft)) may there be expressed in all orders of perturbation theory as one universal 
function with shifted arguments [|39|]. 

It is convenient to represent the change of FF basis order by order for k > 0 as [ |3(~)| | 7 : 

D + {w, jUo) = (1 - aco)D s (o),i4) - E m D g ((o,^), D~((o,/j.q) = a a,D,(m,^)+ e ffl D g (m,^). 

(2.29) 

This implies for the components of the timelike-splitting-function matrix that 

P-l(<o) = acaP^^a) + SaPS {a) + ficoP^ (a) + (1 - a a )P { gg ] (co), 


,(*) 


p?l((0) = p|T(«)- 


(ffl) = P--((0)~ ^P^(0}) + ^-^pff(0}^ , P^l(o})P^(co)+P^(co)-P^l(o}), 


.(*) , 


t(D 


= P£\a)~ (^(ffl)-^(ffl) ) , (2.30) 


>(*) 


Otffl D (k). 


where a 0) , (i 0 > and e w are given in Eq. ( |2. 1 8| ). 

Note, howerver, that the approach ( p .2 8[ ) is not so conveninet in FF case, because we would like 
to keep the diagonal part of P((o,a s ) matrix without an expansion on a s . So. below our approach 
to solve Eq. ( |2.25| ) differs from the usual one (see |[30|l) We write the solution expanding about the 
diagonal part of the all-order timelike-splitting-function matrix in the plus-minus basis, instead of 
its LO contribution. For this purpose, we rewrite Eq. (2.28) in the following way: 


P{(0,a s ) = 


P ++ (m,a s ) 0 
0 P _( co,a s ) 




0 


p+2(w) o 


**-+(«) 


+ 


0 0{a) 

G(a\) 0 


. (2.31) 


In general, the solution to Eq. (2.25|) in the plus-minus basis can be formally written as 


D (P 2 ) = T P- { ex P 


fP 2 Jr ,2 








(2.32) 


where 7b 2 denotes the path ordering with respect to fi 2 and 


D = 



(2.33) 


7 The difference in the diagonalization to compare with the spacelike case considered above is following: 
—2 Pgq and, thus, ft® <->• e ffl . 
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As anticipated, we make the following ansatz to expand about the diagonal part of the timelike- 
splitting-function matrix in the plus-minus basis: 


(exp/ r f 2 p (P 2 )\ -Z '(Miexp 

r ^,> d (p 2 )] 

{ J »Z P~ J 

[M p 


Z(tt 


■Oh 


where 


P D (co) = 


P ++ (co) 0 
0 P—(co) 


(2.34) 


(2.35) 


is the diagonal part of Eq. ( |2.31| ) and Z is a matrix in the plus-minus basis which has a perturbative 
expansion of the form 

Z(jU 2 ) = 1 + a s (fi 2 )Z^ + &(a 2 ). (2.36) 


Changing integration variable in Eq. ( 2.34[ ), we obtain 


T as { exp 


/ 


1 Pm 2 ) da , 


(Mo 2 ) £(«*) 


P{a s )\=Z 1 (fli(jU 2 ))exp 


f a PM~) da 


Wo 2 ) P(a* 


4rP U (a s ) 


Z(a s (p o 2 )). (2.37) 


Substituting then Eq. ( p.36 ) into Eq. (2.37), differentiating it with respect to a s , and keeping only 
the first term in the a s expansion, we obtain the following condition for the Z ( 11 matrix (see Section 
8.1 for the similar procedure in the spacelike case): 


Z (1) + 


Solving it, we find: 


p(0)D 

"AT’ 


p(l)OD 

IT’ 


^ OD (co) = 


P+l(co) 0 




(2.38) 


Z±i(ffl) = 0 , z^(m) = 




Po + Pm<o)-pW(co) 


(2.39) 


At this point, an important comment is in order. In the conventional approach to solve Eq.(2.25), 
one expands about the diagonal LO matrix given in Eq. ( 2.28 ), while here we expand about the all¬ 
order diagonal part of the matrix given in Eq. ( 2.31| ). The motivation for us to do this arises from 
the fact that the functional dependence of P±± ( 0 ). a s ) on a s is different after resummation. 

Now reverting the change of basis specified in Eq. (2.29), we find the gluon and quark-singlet 
fragmentation functions to be given by 


Dg(co,p 1 ) = -^D+( ( o 1 p 2 ) + 


1 — 


D (ft), At 2 ), D cl (co,p 2 ) =D + ((0 ,p 2 )+D (ft), At 2 


(2.40) 

As expected, this suggests to write the gluon and quark-singlet fragmentation functions in the 
following way: 


Da{co,p 2 ) = D+ (w,At 2 ) +D a (to,At 2 


a = g,q, 


(2.41) 


where D+ (ft), At 2 ) evolves like a plus component and D a (o).ji 2 ) like a minus component. 
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We now explicitly compute the functions D±{(0,ji 2 ) appearing in Eq. (2.41). To this end, we 
first substitute Eq. ( 2.34 ) into Eq. ( |2.32 ). Using Eqs. ( 2.35| ) and ( |2.39| ), we then obtain 


D+(w,jU 2 ) = D + (co,ii$)t + (<o,n 2 ,ii$)-a s (n 2 )Z { _l ] + {(o)D (co,^)f_(to,^ 2 ,^), 
D~(co,id 2 ) = D-{(o,^)f_{(o,ix 2 ,^)-a s {^ 2 )Z { ^_{(o)D + {(o,^)f + {(o,^ 2 ,^), (2-42) 


where 


and 


D ± (m,rf))=D ± {m,v 2 )+a s (ii 2 )Z^l(cQ)DT((Q,ii 2 ), 

r fl »(M 2 ) dao 


T±{(o,ix 2 ,^) = ex p 


/ 


(Mo 2 ) Pfa ) 


P±±{(0,a s ) 


(2.43) 


(2.44) 


has a RG-type exponential form. Finally, inserting Eq. Q2.42D into Eq. ( |2.4Q[ ), we find by compari- 

(2.45) 


son with Eq. (2.41) that 


Dt((O,P Z )=Dt{(0,l4)T ± (a),ld 2 ,lj£)H±(a),ii z ), 


where 


(®,Mo 2 ) = -^D+(w,/i 0 2 ), D g = -——-D q (<o,ng), 

to) &CQ 

DU(o,l4)=D + {o),i4), D-{(o,ixl) =D~{(o,ix o 2 ), 


and H^{(o,jX 2 ) are perturbative functions given by 


H±{(0,ix 2 ) = 1 -fl s (ti 2 )zi' i » + % 2 ). 


(2.46) 


(2.47) 


At £?(a s ), we have 


z£l, g {co) = -z2l((o)( 1 ““ 


, Z ( il q {(0)=z2l{(0), (2.48) 

\ Wco J 


where Z±^{oo) is given by Eq. ( [2.39| ). 


3. Generalized DAS approach 


The flat initial condition ( |1.1[ ) corresponds to the case when parton density tend to some 
stant value at x 


0 and at some initial value Q 2 


o- 


The main ingredients of the results [15, 


are: 


con- 


A. Both, the gluon and quark singlet densities are presented in terms of two components (” + ” 
and ” — ”) which are obtained from the analytic 0 2 -dcpcndcnt expressions of the cotTesponding 
(” + ” and ” — ”) PDF moments. 

B. The twist-two part of the ” —” component is constant at small x at any values of Q 2 , 
whereas the one of the ” + ” component grows at Q 2 > Ql as 


r^j £ 


0NLO 


) Cnlo = 2 J [|</+| snlo - (d++ + |d+| b\) pnlo] In , 


Pnlo 


Qnlo 

21n(l/jc) ’ 


(3.1) 
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where a and p are the generalized Ball-Forte variables, 


snlo = In 


\Ql) 


flNL° (e 2 


,, — „NLO/ n 2\ „NL0/ n 2\ a — 

, Pnlo — a s (<J 0 )-a s ), d+ — - 


12 

4 


d++ — 


i “++ 


412 

W 


and /3o and j3i are given in Eq. (2.5) 

3.1 Parton distributions and the structure function F 2 

The results for parton densities and F 2 are following: 

• The structure function F 2 has the form: 

Fl°(x,Q 2 ) = efq„ LO (x,Q 2 ), /a,L 0 (x,e 2 ) = faWW+faWQ 2 ) 

at the LO approximation, where e is the average charge square (|2. 1 1), and 


(3.2) 


(3.3) 


F? lo (x,Q 2 ) = e (/^nloCt Q 2 ) + lfa" LO (Q 2 )f g m.o(x, Q 2 

fa,NLo( x ,Q ) = /oInLO^’G )+/ ai NLo( X )2 ) 

at the NLO approximation. 

The small-* asymptotic results for the LO parton densities f± LO are 

4lo4<2 2 ) = ^A g + ^A^I 0 (a LO ) e~ d+SLO + <9(p LO ), 
Z^WaC? 2 ) = ^ (a 8 + -xAq]p L oh((J L o) e^ +SLO + O(plo), 


where ' 


(3.4) 


(3.5) 

(3.6) 


4lo(*,2 2 ) = ~AA q e- d - s '° + 0(x), /- LO (*,(2 2 ) = A q e~ d - S ^ + 0(x), (3.7) 


-7 _ 20 / t - l6 f 

d-\- — 1 +_, d— — 


(3.8) 


27/3 0 ’ 27 A) 

are the regular parts of the anomalous dimensions d + {n) and d- («), respectively, in the limit 
n —» l 9 . Here n is the variable in Mellin space. The functions I v (v = 0,1) are related to the 
modified Bessel function 7 V and to the Bessel function J v by: 


A(c) = 


hip). 


if 5 > 0 


i v J v (io), ; 2 = — 1, if s < 0 
At the LO, the variables Clo and Plo are given by Eq. ( fO[ ) when p = 0, i.e. 


(3.9) 


Clo — 2w |d+| 5 Loln ( - ) , Plo— 


Clo 


21n(l/x) ’ 

and the variable s L o i s given by Eq. (|3.2|) with a\:°(Q 2 ) as in Eq. (p~T|). 


(3.10) 


8 The dependence on the colour factors C/ = 3, Cf = 4/3 in Eqs. ( |3.2| ), ( [3,8| ) and ( |3. 1 7| ) can be found in [fit]]. 

9 We denote the singular and regular parts of a given quantity k(n) in the limit n-> 1 by k/{n — 1) and k, respectively. 


11 














Parton distributions 


Anatoly Kotikov 


The small-* asymptotic results for the NLO parton densities f+ are 

4 + NLoM 2 ) = A+ NL o(e 2 ,<2o)/o(^NLo)^ +iNLO ' S+PNLO + 0(p NLO), 


/^NLO^jG 2 ) — A+ N LO 


1 — ^+_a^ L0 (G")) Pnlo^i(cnlo) 


+200^° ( Q~)Io( Onlo) 


4nL 0(^2 2 ) =A- NLO (G 2 ,e2)^-( 1 >N L 0-^N L 0 + 

4nlo(^ 2 2 ) =A- NL0 ^- (1) ^o-°-^o + 0(x), 


where (b\ =j8i/j 8 0 ) 
and similar for D + and D + , 


D± = d±± — d±b\ 




/ \ „NLO (n 8 ^/ „NLO/ 


+ 9 ( 1 + ( 3 + ^)< LU (Qo) ~ - gf 4 LU (<2) ) 


'g-NI-O 


(G6Go) =A,-A+ NLO (<26,<2-). 


(3-11) 


e -d + (i) iNLO -D + PNLo + 0 (^ 0 ), (3.12) 


(3.13) 

(3.14) 

(3.15) 


(3.16) 


The coupling constant a s (Q 2 ) is introduced in Eq. (|2.2|). The variables d + , d ++ d + and d 

are diven in Eqs. ( |3.2[ ) and ( |3.8| ), respectively. The nonzero variables d ++ , d _and d" 

(,a = q,g) have the form 


d++ ~ 4l 36C3+33?2_ 'ir + ¥ 


68 13/ 

9 243 


—,g _ 80/ -jg _ 0 / 

’ 61 ’ - +__ 27’ 


sK 26 " 36 ^*' 

with 6 and 6 are Eller functions. 


’ „ 23 13/' 

46--— 

b 18 243 


, 4_ =23-126- 


13/ 

81 


(3.17) 


4. Resummation in FFs 


As already mentioned in Introduction, reliable computations of AJMs require resummed an¬ 
alytic expressions for the splitting functions because one has to evaluate the first Mellin moment 
(corresponding to co = 0), which is a divergent quantity in the fixed-order perturbative approach. 
As is well known, resummation overcomes this problem, as demonstrated in the pioneering works 


by Mueller [ ]18Q and others [[4(J |. 

In particular, as we shall see in previous subsection, resummed expressions for the first Mellin 


moments of the timelike splitting functions in the plus-minus basis appearing in Eq. (2.28) are 
required in our approach. Up to the NNLL level in the MS scheme, these may be extracted from 
the available literature flT8|, p~9|. |2C|. 21 ] in closed analytic form using the relations in Eq. ([2 .3( ). 
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For future considerations, we remind the reader of an assumpion already made in Ref. [20] 
according to which the splitting functions P__(<o) and are supposed to be free of sin¬ 

gularities in the limit co —>• 0. In fact, this is expected to be true to all orders. This is certainly 
true at the LL and NLL levels for the timelike splitting functions, as was verified in |[2()||. This is 
also true at the NNLL level, as may be explicitly checked by inserting the results of Ref. [13] in 
Eq. ( 2.30| ). Moreover, this is true through NLO in the spacelike case [(T5J] and holds for the LO and 
NLO singularities [ pdj . [42[ . |39| ] to all orders in the framework of the BFKL dynamics ||sj], a fact that 


was exploited in various approaches (see, e.g., Refs. [43] and references cited therein). We also 
note that the timelike splitting functions share a number of simple properties with their spacelike 
counterparts. In particular, the LO splitting functions are the same, and the diagonal splitting func¬ 
tions grow like In co for co —> °° at all orders. This suggests the conjecture that the double-logarithm 
resummation in the timelike case and the BFKL resummation in the spacelike case are only related 
via the plus components. The minus components are devoid of singularities as co —>• 0 and thus are 
not resummed. Now that this is known to be true for the first three orders of resummation, one has 
reason to expect this to remain true for all orders. 

Using the relationships between the components of the splitting functions in the two bases 

given in Eq. ( |2.30 ), we find that the absence of singularities for co = 0 in P _ ((0,a s ) and P+-(<o,a s ) 

implies that the singular terms are related as 


P^ g (co,a s ) = ——^-Ppl ng (co,a s ), P™ g (( 0 ,a s ) = P™ g (co,a s ), 


OC a 


(4.1) 


where, through the NLL level, 


10 


£(o 


Ca 

C f 


i-“ 

6 


1+2 Z>_ 4 q^\ 

C, Cl 


) 


+ < 


co- 


(4.2) 


An explicit check of the applicability of the relationships in Eqs. ( PI ) for Pij{co,a s ) with ij = g,g 
themselves is performed in the Appendix of Ref. [23]. Of course, the relationships in Eqs. (0) 
may be used to fix the singular terms of the off-diagonal timelike splitting functions P qg (co,a s ) 
and P gq (co,a s ) using known results for the diagonal timelike splitting functions P qq (co,a s ) and 
P gg ((0,a s ). Since Refs. [19, 38] became available during the preparation of Ref. [20], the relations 
in Eqs. ( p~T| ) provided an important independent check rather than a prediction. 

We take here the opportunity to point out that Eqs. (2.45) and ( 2.46| ) together with Eq. ( |4.2| ) 
support the motivations for the numerical effective approach that we used in Ref. [28, [231] to study 


the gluon-to-quark AJM ratio. In fact, according to the findings of Ref. [28, 23], substituting 
co = ft) e rr> where 

= 2^2 Caci s , (4.3) 


into Eq. ( fk2| ) exactly reproduces the result for the average gluon-to-quark jet multiplicity ratio 
r(Q 2 ) obtained in Ref. [|44||. In the next section, we shall obtain improved analytic formulae for the 
ratio r(Q 2 ) and also for the average gluon and quark jet multiplicities. 

1(l To have a possibility to compare different approximations, it is convenient to keep the general forms of the colour 
factors Ca — 3,Cf — 4/3 in the present and the next sections. 
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Here we would also like to note that, at first sight, the substitution to = to e g- should induce a 
Q 2 dependence to the diagonalization matrix. This is not the case, however, because to double- 
logarithmic accuracy the Q 2 dependence of a s (Q 2 ) can be neglected, so that the factor a m /e m does 
not recieve any Q 2 dependence upon the substitution 0 ) = to e ff. This supports the possibility to 
use this substitution in our analysis and gives an explanation of the good agreement with other 
approaches, e.g. that of Ref. Q]. Nevertheless, this substitution only carries a phenomenological 
meaning. It should only be done in the factor a 0) /e 0J , but not in the RG exponents of Eq. ( 2.44 ), 
where it would lead to a double-counting problem. In fact, the dangerous terms are already re¬ 
summed in Eq. (2.44). 

In order to be able to obtain the AJMs, we have to first evaluate the first Mellin momoments of 
the timelike splitting functions in the plus-minus basis. According to Eq. ( |2.30| ) together with the 
results given in Refs. [18, 21 ], we have 


d nnll 

++ 


(® = 0) = 70(1 — ^70 + K 2 yfi), 


(4.4) 


where 


7o = ^++( w = 0) = s/2C A a s , Ki = 


1 

12 


11+4 


II 

Ca 


1 - 


2 C F 

cT 


Kl 288 


1193-576^-56^ 

C A 




+ 16 # 

+4 


,C f 


.C 


1 +4—— 12—E 

Ca Ci 


and 


where 


-jNNLL 


(ft) = 0) = — 


Cf 

Ca 


pNNLL 
qg 


(o = 0), 


P" nll (cd = 0) = —7>+ --T f 


16. 

T 


2 

3 : 


17-4 


II 

Ca 


1- 


2 C F 

cT 


(2 C A a]) 


1/2 


For the P+- component, we obtain 


pNNLL 

4— 


(to = 0) = &{a 2 s ). 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


Finally, as for the P _component, we note that its LO expression produces a finite, nonvanishing 

term for to = 0 that is of the same order in a s as the NLL-resummed results in Eq. (fi~4|), which 
leads us to use the following expression for the P _component: 


pNNLL^_ r\\ _ 8 TfCf s?r„2\ 

— (CO-0)-- 


(4.10) 


at NNLL accuracy. 


We can now perform the integration in Eq. (2. 


t NNLL/' / q2\ 
nNNLL /a rP- /->2\ 7 ± (G ) 


J (0,Q~,Qo) = 


through the NNLL level, which yields 

NNLL / rP \ t i NLL/ / a2' / , ~ts\d. 


j’NNLL ( rP 


r NNLL 

' + 


( Q') = exp 


(Go 

4Ca 

W«2 




(Q Z ) = T^ l (Q a ) = (a s {Q 2 )) 

[1 + (Zq -2C A K 2 )a s (Q 2 )\ \ (, a s (Q 2 )) d+ 


where 


/ & 
bl = W 


d _ = 


8 T f C f 
3Ca(3 0 ’ 


d + = 


2C a K\ 

Po 


(4.11) 

(4.12) 

(4.13) 
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5. Multiplicities 


According to Eqs. ( 2.44 ) and ( 2.45 ), the components are not involved in the AJM Q 2 
evolution, which is performed at 0 ) = 0 using the resummed expressions for the plus and minus 
components given in Eq. ( |4.4[ ) and ( 4.1 Op , respectively. We are now ready to define the average 
gluon and quark jet multiplicities in our formalism, namely 


(nh(Q 2 )) a = D a (0,Q 2 ) = D^{0,Q 2 )+D~(0,Q 2 ), (a = g,q). 

On the other hand, from Eqs. ( |2.45| ) and ( [2.46| ), it follows that 


(5.1) 


, D+(0,Q 

r +(Q 2 ) = 


« 0 ) H g ( 0 ) -Q 2 
D q (0, Q 2 ) w'Ab Hq((o,Q 2 


= - fim — y r (<2 2 ) = ° 8 (Q ’ g ) = lim 1 aa Hg ^ ,Q ^ 

r + t,,n2V -W) Dq{ 0)G 2 ) H q (co,Q 2 ) 


(5.2) 

Using these definitions and again Eq. ( |2.45| ), we may write general expressions for the gluon and 
quark AJMs: 


(n h (Q-))g = D+(0, <2^)r; es (0, Q\Qo)Hg(QjQ) +D~{ 0, Qb)r-(Q-)T™(0, Q\Q Z 0 )H~( 0, Q 1 


{n h (.Q 2 ))s = 


Dg{0,Ql) 

r+m 


r; es (0, Q 2 ,Qq)H+(0,Q 2 ) + D-(0, <25)r res (0, Q 2 ,Ql)H q (0 ,Q 


At the LO in a s , the coefficients of the RG exponents are given by 

r+ (Q 2 ) = ^, r_(<2 2 )=0, // i ± (0,e 2 ) = l, D±(0,Q 2 )=D±(0,Q 2 ). 

Cf 


(5.3) 


(5.4) 


It would, of course, be desirable to include higher-order corrections in Eqs. ( p.4[ ). However, 
this is highly nontrivial because the general perturbative structures of the functions (o). q 2 ) 
and Z± TtCI ((o,a s ), which would allow us to resum those higher-order corrections, are presently 
unknown. Fortunatly, some approximations can be made. On the one hand, it is well-known that 
the plus components by themselves represent the dominant contributions to both the gluon and 
quark AJMs (see, e.g.. Ref. [ j45| | for the gluon case and Ref. [46] for the quark case). On the other 
hand, Eq. (5.2) tells us that D~(0,Q 2 ) is suppressed with respect to D c/ (0. Q 1 ) because a (0 ~ I + 


These two observations suggest that keeping r .AQ 2 ) = 0 also beyond LO should represent 
a good approximation. Nevertheless, we shall explain below how to obtain the first nonvanishing 
contribution to t\ (Q 2 ). Furthermore, we notice that higher-order corrections to H+((). Q 2 ) and 
Df (0, Qq) just represent redefinitions of (0, <2o) by constant factors apart from running-coupling 
effects. Therefore, we assume that these corrections can be neglected. 


Note that the resummation of the ±± components was performed similarly to Eq. ( |2.44| ) for 
the case of parton distribution functions in Ref. & Such resummations are very important be¬ 
cause they reduce the Q 2 dependences of the considered results at fixed order in perturbation theory 
by properly taking into account terms that are potentially large in the limit to 


0 [47, 


. We 

anticipate similar properties in the considered case, too, which is in line with our approximations. 
Some additional support for this may be obtained from jV = 4 super Yang-Mills theory, where the 
diagonalization can be performed exactly in any order of perturbation theory because the coupling 
constant and the corresponding martices for the diagonalization do not depended on Q 2 . Conse¬ 


quently, there are no Z±L a ((Q) terms, and only i±±(®) terms contribute to the integrand of the RG 
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exponent. Looking at the r.h.s. of Eqs. (2.43) and (2.47), we indeed observe that the corrections of 
&{a s ) would cancel each other if the coupling constant were scale independent. 

We now discuss higher-order corrections to r + (Q 2 ). As already mentioned above, we intro¬ 


duced in Ref. [28] an effective approach to perform the resummation of the first Mellin moment of 
the plus component of the anomalous dimension. In that approach, resummation is performed by 
taking the fixed-order plus component and substituting co = <y e ff, where a) e ff is given in Eq. (fi~3|). 
We now show that this approach is exact to G{^fa s ). We indeed recover Eq. ( fP] ) by substituting 
0 ) = ft) e ff in the leading singular term of the LO splitting function P++(co,a s ), 


= 


4 C A a s 


(0 




(5.5) 


We may then also substitute co = ft) e ff in Eq. (5.2) before taking the limit in 0 ) = 0. Using also 
Eq. (0), we thus find 


MQ 2 ) = 


Ca 

C f 


1 - 




1 + 2 


T f C f T f \ 


C a 


C 2 / 


+ C?(a s ), 


(5.6) 


which coincides with the result obtained by Mueller in Ref. [44]. For this reason and because, in 


Ref. [}49[] , the gluon and quark AJMs evolve with only one RG exponent, we inteprete the result in 
Eq. (5) of Ref. [^5]j as higher-order corrections to Eq. ([+?]). Complete analytic expressions for all 


the coefficients of the expansion through G(a s ) may be found in Appendix 1 of Ref. [[25]]. This 
interpretation is also explicitely confirmed in Chapter 7 of Ref. [ |50| ] through (Ha s ). 

Since we showed that our approach reproduces exact analytic results at we may safely 

apply it to predict the first non-vanishing correction to r_(<2 2 ) defined in Eq. (|+2|), which yields 


(< 2 2 ) = -^ 


47> 2a s {Q 2 ) 


C A 


+ &(a s ). 


(5.7) 


However, contributions beyond ff{yfoT s ) obtained in this way cannot be trusted, and further inves¬ 
tigation is required. Therefore, we refrain from considering such contributions here. 

For the reader’s convenience, we list here expressions with numerical coefficients for r + (Q 2 ) 


through 


3/2 \ 


and for r_ (Q 2 ) through in QCD with nj = 5: 


r + (Q 2 ) = 2.25-2.m49Ja s {Q 2 )-n.54a s (Q 2 ) + lQM62a s /2 (Q 2 ) + 0(a 


r~(Q 2 ) = —2.72166 ^ a s (Q 2 ) + ff(a s ). 


(5.8) 

(5.9) 


We denote the approximation in which Eqs. ( [1. 1 l| )-( p~l2| ) and (5.4) are used as LO + NNLL, 
the improved approximation in which the expression for r + (Q 2 ) in Eq. ( +4) is replaced by Eq. (+8), 
i.e. Eq. (5) in Ref. [[25)], as N 3 LO a p prox + NNLL, and our best approximation in which, on top of 
that, the expression for r~(Q 2 ) in Eq. (5.4) is replaced by Eq. 


as 


N 3 LO 


approx 


+ NLO + NNLL. 

We shall see in the next Section, where we compare with the experimental data and extract the 
strong-coupling constant, that the latter two approximations are actually very good and that the last 
one yields the best results, as expected. 
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In all the approximations considered here, we may summarize our main theoretical results for 
the gluon and quark AJMs in the following way: 


(n h (Q 2 )) g = n i (£?o) ?+ es (0, 2 2 , Go) + «2 (<2o) r -(Q 2 ) TL es (0, £> 2 , Go), 


(n h (Q 2 )) s = m (So) r " S + "2 ( 65 ) f - es (0, Q 2 , Q 2 


r+(Q 2 


0 T 


where 


.. .. ,„2s D MQo)-r-(Q 2 o)D s (0,Q 2 o) .. ^ 2 , r + (G 2 ) Di(0 , ! 22 ) _ Dg( o, ! 22 ) 

'filGoJ = r +(Go) - /^ 2 \ — -> n 2 {Qo) = 


r + {Ql)-r.{Ql) 

The gluon-to-quark AJM ratio may thus be written as 


r + {Ql)-r.{Ql) 


(5.10) 


(5.11) 


Ke 2 ) = m^ = r + (e 2 ) 


(nh(Q 2 )) s 


2 m 


where 


l + r^Q 2 )R(Q 2 )T™(0,Q 2 ,Q 2 )/T™(0,Q 2 ,Q 2 
ll + r + (Q 2 )R(Q 2 )T™(0,Q 2 ,Q 2 )/f?s(0,Q 2 ,Q 2 )_ 

ml) =" 


>mi) 

«l(Gn) 


(5.12) 


(5.13) 


It follows from the definition of 7± es (0, Q 2 ,Qq) in Eq. ( |4. 1 1| ) and from Eq. (5.11) that, for Q 2 = Qfc, 
Eqs. (5.10) and (5.12) become 

(n h (Q 2 o)) g = D g (0,Qo), (n h (Q 2 0 )) q =D s (0,Q 2 ), r(Q §) = 


p g (ml) 

D s ( 0,Q 2 ) 


(5.14) 


These represent the initial conditions for the Q 2 evolution at an arbitrary initial scale <2o- I n fact, 
Eq. ( [5. 1 (] ) is independ of Q 5 , as may be observed by noticing that 

tr(0,Q 2 ,Q 2 0 ) = fr(0,Q 2 ,Qi)fr(0,Qj,Ql), (5.15) 

for an arbitrary scale <2i (see also Ref. [ pTj | for a detailed discussion of this point). 

In the approximations with r_(<2 2 ) = 0 [22], i.e. the LO + NNLL and N 3 LO approx + NNLL 
ones, our general results in Eqs. ( 5. 10| ), and ( 5.12 ) collapse to 


(n h (Q 2 )) g = D g (0,Q 2 )T™(0,Q 2 ,Ql), 

21 ) 


(n h (Q 2 )) s =D g (0,Qo) ±+ 

m 2 ) = 


fr(0,Q 2 ,Q 2 ^ 


r + (Q 2 ) 

r +(Q 2 ) 


D s (0,Ql)~ 


DMQl 

r + (Ql) 


TL es (0,Q 2 ,Q 2 ), 


1 I r+(Q 2 ) 

1 + MG6) 


D s {0,Ql)r+(Ql) 

D g (0,Q 2 ) 


T™(0,Q 2 ,Ql) 

TF(0,Q 2 ,Q 2 0 ) 


(5.16) 


The NNLL-resummed expressions for the gluon and quark AJMs given by Eq. ( 5.1 0| ) only 
depend on two nonperturbative constants, namely D g (0,Qfi) and D s (0,Qq). These allow for a 
simple physical interpretation. In fact, according to Eq. ( |5. 14 ), they are the average gluon and quark 
jet multiplicities at the arbitrary scale <2o- We should also mention that identifying the quantity 
r \(Q 2 ) with the one computed in Ref. [p5[], we assume the scheme dependence to be negligible. 
This should be justified because of the scheme independence through NLL established in Ref. [EH]. 
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Figure 1: F 2 (x,Q 2 ) as a function of x for different Q 1 bins. The experimental points are from HI ||T] (open 
points) and ZEUS ||]| (solid points) at Q 2 > 1.5 GeV 2 . The solid curve represents the NLO fit. The dashed 
curve (hardly distinguishable from the solid one) represents the LO fit. 


We note that the Q 2 dependence of our results is always generated via a s {Q 2 ) according to 
Eq. (|2.4|). This allows us to express Eq. ( [4.1 l[ ) entirely in terms of oc s (Q 2 ). In fact, substituting the 
QCD values for the color factors and choosing rif = 5 in the formulae given in Refs. [22, 23], we 
may write at NNLL 


\Q\Ql) = 


T± S (Q 2 ) 

Ti es (Q 2 0 )' 


T™{Q 2 ) = of 1 (Q 2 ), r[ es (<2 2 ) = exp 


di + dj, cx s (Q 2 ) 

\M(2 2 ) 


af 4 (<2 2 ), 


(5.17) 


where 


d { =0.38647, d 2 =2.65187, d 3 = -3.87674, d 4 = 0.97771. (5.18) 
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Table 1: The result of the LO and NLO fits to HI and ZEUS data for different low Q 2 cuts. In the fits / is 
fixed to 4 flavors. 




4? 

Qi [GeV 2 ] 

X 2 /n.o.p. 

Q 2 > 1.5 GeV 2 
LO 

LO&an. 

LO&fr. 

0.784±.016 

0.932±.017 

1.022±.018 

0.801±.019 

0.707±.020 

0.650±.020 

0.304±.003 

0.339±.003 

0.356±.003 

754/609 

632/609 

547/609 

NLO 

NLO&an. 

NLO&fr. 

-0.200±.011 

0.310±.013 

0.180±.012 

0.903±.021 

0.640±.022 

0.780±.022 

0.495±.006 

0.702±.008 

0.661±.007 

798/609 

655/609 

669/609 

Q 2 > 0.5GeV 2 
LO 

LO&an. 

LO&fr. 

0.641±.010 

0.846±.010 

1.127±.011 

0.937±.012 

0.771±.013 

0.534±.015 

0.295±.003 

0.328±.003 

0.358±.003 

1090/662 

803/662 

679/662 

NLO 

NLO&an. 

NLO&fr. 

-0.192±.006 

0.281±.008 

0.205±.007 

1.087±.012 

0.634±.016 

0.650±.016 

0.478±.006 

0.680±.007 

0.589±.006 

1229/662 

633/662 

670/662 


6. Comparison with experimental data 


Here we compare our formulae with experimental data for DIS SF /^(x, Q 2 ) and for the AJMs. 
In the DIS case, we linrite ourselves by consideration only the SF F 2 (x, Q 2 ). The comparison of 
the generalized DAS approach predictions with the data for the slope r) In File) 1 n(1 /x) [|4]| and for 
the heavy parts of F% [H] can be found in Refs. [52, Iff]] and [53], respectively (see also the review 


p). An estimation of the cross-sections of very high-energy neutrino and nucleon scattering has 
been found in [|55|]. 

6.1 DIS SF F 2 

Using the results of section 3 we have analyzed HERA data for Fo at small x from the HI and 
ZEUS Collaborations [|I|, [2], |j]. 

In order to keep the analysis as simple as possible, we fix / = 4 and a s (M |) = 0.1166 (i.e., 
A^ 4 ^ = 284 MeV) in agreement with the more recent ZEUS results |||]. 

As it is possible to see in Fig. 1 (see also [@-| [T6|] ), the twist-two approximation is reasonable 
at Q 2 > 2 GeV 2 . At smaller Q 2 , some modification of the approximation should be considered. 


In Ref. [ |1 6p we have added the higher twist corrections. For renormalon model of higher twists, 
we have found a good agreement with experimental data at essentially lower Q 2 values: Q 2 >0.5 
GeV 2 (see Figs. 4 and 5 in Ref. [jlfj), but we have added 4 additional parameters: amplitudes of 
twist-4 and twist-6 corrections to quark and gluon densities. 

Moreover, the results of fits in [|T^] have an important property: they are very similar in LO and 
NLO approximations of perturbation theory. The similarity is related to the fact that the small-x 
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Figure 2: x dependence of I-'iix- Q 2 ) in bins of Q 2 . The experimental data from HI (open points) and ZEUS 
(solid points) are compared with the NLO fits for Q 2 > 0.5 GeV 2 implemented with the canonical (solid 
lines), frozen (dot-dashed lines), and analytic (dashed lines) versions of the strong-coupling constant. For 
comparison, also the results obtained in Ref. [[l6]] through a fit based on the renormalon model of higher-twist 
terms are shown (dotted lines). 


asymptotics of the NLO corrections are usually large and negative (see, for example, a s -corrections 
pT| , [42]] to BFKL kernel [^] 1 ’). Then, the LO form ~ oc s (Q 2 ) for some observable and the NLO 
one ~ a s (Q 2 )(\ — KcxjQ 2 )) with a large value of K are similar, because Anlo Alo 12 and, thus, 
a s (Q 2 ) at LO is considerably smaller then a s ( Q 2 ) at NLO for HERA Q 2 values. 

In other words, performing some resummation procedure (such as Grunberg’s effective-charge 
method [|6|]), one can see that the results up to NLO approximation may be represented as ~ 
cts(<2 2 ff)’ where Q 2 IT Q 2 . Indeed, from different studies [|57|, 58, 59], it is well known that at 


11 It seems that it is a property of any processes in which gluons, but not quarks play a basic role. 

12 The equality of a s (M |) at LO and NLO approximations, where Mz is the Z-boson mass, relates Anlo an d Alo: 
A^q = 284 MeV (as in [jj]) corresponds to Alo = 112 MeV (see ||I^] ). 
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small-x values the effective argument of the coupling constant is higher then Q 2 . As it was shown 
in [^], the usage of the effective scale in the generalized DAS approach improves the agreement 
with data for SF F 2 (x, Q 2 ). 

Here, to improve the agreement at small Q 2 values without additional parameters, we mod¬ 
ify the QCD coupling constant. We consider two modifications, which effectively increase the 
argument of the coupling constant at small Q 2 values (in agreement with [57[ 58, 59]). 

In one case, which is more phenomenological, we introduce freezing of the coupling constant 
by changing its argument Q 2 —>• Q 2 +M 2 , where M p is the p-meson mass (see [|5l|]). Thus, in the 
formulae of the Section 2 we should do the following replacement: 


a s (Q)^ ail (Q 2 )=a s (Q 2 +M-) 


( 6 . 1 ) 


The second possibility incorporates the Shirkov-Solovtsov idea [p2|]-[p5|] about analyticity of 
the coupling constant that leads to the additional its power dependence. Then, in the formulae of 
the previous section the coupling constant a s (Q 2 ) should be replaced as follows: 


a^(Q z )=a s (Q 2 )~ 


1 


^■LO 


A) Q 2 - At 


( 6 . 2 ) 


LO 


at the LO approximation and 


flan(G 2 ) = a s {Q 2 )- 


1 


A 2 


■ + ..., 


(6.3) 


2A) Q 2 - A 2 

at the NLO approximation, where the symbol ... stands for terms which have negligible contribu¬ 


tions at Q > 1 GeV @] 13 . 

Figure 2 and Table 1 show a strong improvement of the agreement with experimental data for 
F 2 (y 2 values decreased almost 2 times!). 


6.1.1 Hl&ZEUS data 

Here we have analyzed the very precise Hl&ZEUS data for F 2 ]^|. As can be seen from Fig. 3 
and Table 2, the twist-two approximation is reasonable for Q 2 >4 GeV 2 . At lower Q 2 we observe 
that the fits in the cases with “frozen” and analytic strong coupling constants are very similar (see 
also [66]) and describe the data in the low Q 2 region significantly better than the standard fit (y 2 
values decreased 2 a 3 times!) Nevertheless, for Q 2 < 1.5 GeV 2 there is still some disagreement 
with the data, which needs to be additionally studied. In particular, the BFKL resummation [pi] 


may be important here [67]. It can be added in the generalized DAS approach according to the 


discussion in Ref. [54], 


6.2 Average multiplicity and experimendal data 


Now we show the results in [|23|] obtained from a global fit to the available experimental data 
of our formulas in Eq. ( 5.10| ) in the LO + NNLL, N 3 LO ap prox + NNLL, and N 3 LO appr0 x + NLO + 


13 Note that in inm more accurate, but essentially mor e cumbersome approximations of a an (Q 2 ) have been 
proposed. We limit ourselves by above simple form ( |6.3| ) and plan to add the other modifications in our future 
investigations. 
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F 2 (x,Q 2 ) 



Figure 3: x dependence of ^(x, Q 1 ) in bins of Q 2 . The combined experimental data from HI and ZEUS 
Collaborations |||] are compared with the NLO fits for Q 2 > 0.5 GeV 2 implemented with the standard (solid 
lines), frozen (dot-dashed lines), and analytic (dashed lines) versions of the strong coupling constant. 

NNLL approximations, so as to extract the nonperturbative constants D g (0, Qq) and D s (0, Q^). 
We have to make a choice for the scale Qq, which, in principle, is arbitrary. In [[23]], we adopted 
Qo = 50 GeV. 

The gluon and quark AJMs extracted from experimental data strongly depend on the choice 
of jet algorithm. We adopt the selection of experimental data from Ref. [^8|] performed in such 
a way that they correspond to compatible jet algorithms. Specifically, these include the gluon 
AJM measurements in Refs. [|68[]-[|72|] and quark ones in Refs. [ |69| , |73| ], which include 27 and 51 
experimental data points, respectively. The results for (nh{Qo)) g and (»/,(<2o))<? at Qo = 50 GeV 
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Table 2: The results of LO and NLO fits to HI & ZEUS data [|]], with various lower cuts on Q 1 : in the fits 
the number of flavors / is fixed to 4. 



Ag 

At 

Q 2 [GeV 2 ] 

X 2 /n.d.f. 

Q 2 > 5 GeV 2 





LO 

0.623T0.055 

1.204i0.093 

0.437i0.022 

1.00 

LO&an. 

0.796T0.059 

1.103i0.095 

0.494i0.024 

0.85 

LO&fr. 

0.782T0.058 

l.U0±0.094 

0.485i0.024 

0.82 

NLO 

-0.252T0.041 

1.335i0.100 

0.700±0.044 

1.05 

NLO&an. 

0.102±0.046 

1.029i0.106 

1.017i0.060 

0.74 

NLO&fr. 

-0.132±0.043 

1.219i0.102 

0.793i0.049 

0.86 

Q 2 > 3.5GeV 2 





LO 

0.542i0.028 

1.089i0.055 

0.369i0.011 

1.73 

LO&an. 

0.758±0.031 

0.962i0.056 

0.433i0.013 

1.32 

LO&fr. 

0.775±0.031 

0.950i0.056 

0.432i0.013 

1.23 

NLO 

-0.310±0.021 

1.246±0.058 

0.556i0.023 

1.82 

NLO&an. 

0.116T0.024 

0.867i0.064 

0.909±0.330 

1.04 

NLO&fr. 

-0.135i0.022 

1.067i0.061 

0.678i0.026 

1.27 

Q 2 > 2.5GeV 2 





LO 

0.526±0.023 

1.049i0.045 

0.352i0.009 

1.87 

LO&an. 

0.761i0.025 

0.919i0.046 

0.422i0.010 

1.38 

LO&fr. 

0.794±0.025 

0.900±0.047 

0.425i0.010 

1.30 

NLO 

-0.322i0.017 

1.212i0.048 

0.517i0.018 

2.00 

NLO&an. 

0.132i0.020 

0.825i0.053 

0.898i0.026 

1.09 

NLO&fr. 

-0.123i0.018 

1.016±0.051 

0.658i0.021 

1.31 

Q 2 > 0.5GeV 2 





LO 

0.366±0.011 

1.052i0.016 

0.295i0.005 

5.74 

LO&an. 

0.665±0.012 

0.804i0.019 

0.356i0.006 

3.13 

LO&fr. 

0.874T0.012 

0.575i0.021 

0.368i0.006 

2.96 

NLO 

-0.443i0.008 

1.260±0.012 

0.387i0.010 

6.62 

NLO&an. 

0.121i0.008 

0.656±0.024 

0.764i0.015 

1.84 

NLO&fr. 

-0.071i0.007 

0.712i0.023 

0.529i0.011 

2.79 



LO + NNLL 

N 3 LO a pp rox ±NNLL 

N 3 LO approx ± NLO ± NNLL 

(MQl))g 
( n h {Ql)) q 

*dof 

24.31 ±0.85 

15.49 ±0.90 

18.09 

24.02 ±0.36 

15.83 ±0.37 

3.71 

24.17 ±0.36 

15.89±0.33 

2.92 


Table 3 : Fit results for (rc/,(2o))g an d ( n h(Qo))q at Qo = 50 GeV with 90% CL errors and minimum values of 
achieved in the LO + NNLL, N 3 LO a pprox + NNLL, and N 3 LOap prox + NLO + NNLL approximations. 


together with the %% o{ values obtained in our LO + NNLL, N 3 LO appr ox ± NNLL, and N 3 LO appr0 x ± 
NLO + NNLL fits are listed in Table |3| The errors correspond to 90% CL as explained above. All 
these fit results are in agreement with the experimental data. Looking at the %^ of values, we observe 
that the qualities of the fits improve as we go to higher orders, as they should. The improvement is 
most dramatic in the step from LO + NNLL to N 3 LO appr ox + NNLL, where the errors on (n/, ( Qo)) g 
and (n h {Ql)) q are more than halved. The improvement in the step from N 3 LO appr ox ± NNLL 
to N 3 LO appr ox + NLO + NNLL, albeit less pronounced, indicates that the inclusion of the first 
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Figure 4: The average gluon (upper curves) and quark (lower curves) jet multiplicities evaluated from Eq. (5.10), 
respectively, in the LO + NNLL (dashed/gray lines) and N 3 LO appl - 0 x + NLO + NNLL (solid/orange lines) approxima¬ 
tions using the corresponding fit results for (u/,((2o)}g an d ( n h{Qo))q from Table 3 are compared with the experimental 
data included in the fits. The experimental and theoretical uncertainties in the N 3 LO approx + NLO + NNLL results are 
indicated by the shaded/orange bands and the bands enclosed between the dot-dashed curves, respectively. 


correction to r_ (Q 2 ) as given in Eq. ( |5.7| ) is favored by the experimental data. We have verified that 
the values of %% o{ are insensitive to the choice of <2o, as they should. Furthermore, the central values 
converge in the sense that the shifts in the step from N J LO approx + NNLL to N 3 LO approx + NLO + 
NNLL are considerably smaller than those in the step from LO + NNLL to N 3 LO approx + NNLL 
and that, at the same time, the central values after each step are contained within error bars before 
that step. In the fits presented so far, the strong-coupling constant was taken to be the central value 
of the world avarage, af\m\) = 0.1184 [[74]]. In the next Section, we shall include a[ 5 \m\) 
among the fit parameters. 


In Fig. [f[ we show as functions of Q the gluon and quark AJMs evaluated from Eq. ( ]5.10| ) 
at LO + NNLL and N 3 LO approx + NLO + NNLL using the corresponding fit results for («/,(<2o))g 
and (nh(Q})))q at <2o = 50 GeV from Table ||] For clarity, we refrain from including in Fig. [| the 
N 3 LO approx + NNLL results, which are very similar to the N 3 LO approx + NLO + NNLL ones already 
presented in Ref. [22]. In the N 3 LO approx + NLO + NNLL case, Fig. [| also displays two error 
bands, namely the experimental one induced by the 90% CL errors on the respective fit parameters 
in Table |3] and the theoretical one, which is evaluated by varying the scale parameter between (9/2 
and 22- 


While our fits rely on individual measurements of the gluon and quark AJMs, the experimental 
literature also reports determinations of their ratio; see Refs. [[27|. ^8, 70, 72], [75|], which essentially 
cover all the available measurements. In order to find out how well our fits describe the latter 
and thus to test the global consistency of the individual measurements, we compare in Fig. || the 
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Figure 5: The average gluon-to-quark jet multiplicity ratio evaluated from Eq. (5.12) in the LO + NNLL (dashed/gray 
lines) and N 3 LO apP rox + NLO + NNLL (solid/orange lines) approximations using the corresponding fit results for 
( n h(Qo)}g an d ( n h(Qo))q from Table 3 are compared with experimental data. The experimental and theoretical un¬ 
certainties in the N 3 LO appro x + NLO + NNLL result are indicated by the shaded/orange bands and the bands enclosed 
between the dot-dashed curves, respectively. The prediction given by Eq. (5.8) | [f5| l is indicated by the continuous/gray 
line. 


experimental data on the gluon-to-quark AJM ratio with our evaluations of Eq. ( 5. 12| ) in the LO + 
NNLL and N 3 LO a pp 10 x + NLO + NNLL approximations using the corresponding fit results from 
Table |j. As in Lig. [| we present in Lig. || also the experimental and theoretical uncertainties in 
the N 3 LO appr ox + NLO + NNLL result. Lor comparison, we include in Lig. |5] also the prediction of 


Ref. [ )25p given by Eq. 

Looking at Lig. |], we observe that the experimental data are very well described by the 
N ’LOappiox + NLO + NNLL result for Q values above 10 GeV, while they somewhat overshoot 
it below. This discrepancy is likely to be due to the fact that, following Ref. [|68|], we excluded the 
older data from Ref. [27] from our fits because they are inconsistent with the experimental data 
sample compiled in Ref. 


The Monte Carlo analysis of Ref. [ ]26[ ] suggests that the average gluon and quark jet mul¬ 
tiplicities should coincide at about Q = 4 GeV. As is evident from Lig. ||, this agrees with our 
N 3 LO approx + NLO + NNLL result reasonably well given the considerable uncertainties in the 
small-g 2 range discussed above. 

As is obvious from Lig. |5| the approximation of r((Q 2 ) by r + (Q 2 ) given in Eq. (5.8) [25] leads 
to a poor approximation of the experimental data, which reach up to Q values of about 50 GeV. It 
is, therefore, interesting to study the high-2 2 asymptotic behavior of the average gluon-to-quark 
jet ratio. This is done in Lig. where the N 3 LO appr ox + NLO + NNLL result including its exper¬ 
imental and theoretical uncertainties is compared with the approximation by Eq. (|5.8|) way up to 
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Ve 7 


Figure 6: High-g extension of Fig. 5. 



N 3 LO a p prox + NNLL 

N 3 LO approx ± NLO ± NNLL 

(MQo)) g 

{nh(Ql))q 

a s (5) (m|) 

*dof 

24.18 ±0.32 

15.86±0.37 

0.1242 ±0.0046 

2.84 

24.22 ±0.33 

15.88 ±0.35 

0.1199 ±0.0044 

2.85 


Table 4: Fit results for (rc/,(<2o))g and (n/i(Go))? at Go = 50 GeV and for aj 5 ^ (m|) with 90% CL errors and minimum 
values of achieved in the N 3 LO approx + NNLL and N^LOapprox + NLO + NNLL approximations. 

Q = 100 TeV. We observe from Fig. [6] that the approximation approaches the N 3 LO appr ox ±NLO± 
NNLL result rather slowly. Both predictions agree within theoretical errors at Q = 100 TeV, which 
is one order of magnitude beyond LHC energies, where they are still about 10% below the asymp¬ 
totic value Ca/Cf = 2.25. 

6.2.1 Determination of strong-coupling constant from average multiplicity 

In the previous Section, we took ots (. m |) to be a fixed input parameter for our fits. Motivated 
by the excellent goodness of our N 3 LO approx + NNLL and N 3 LO appr0 x + NLO + NNLL fits, we 
now include it among the fit parameters, the more so as the fits should be sufficiently sensitive 
to it in view of the wide Q 2 range populated by the experimental data fitted to. We fit to the 
same experimental data as before and again put <2o = 50 GeV. The fit results are summarized in 
Table [| We observe from Table |] that the results of the N 3 LO appr0 x ± NNLL [ [5~i~| | and N 3 LO appr ox + 
NLO + NNLL fits for («/,(<2o))g an d (nh{Qo))q are mutually consistent. They are also consistent 
with the respective fit results in Table [|. As expected, the values of ;zj o| - are reduced by relasing 
aP (m|) in the fits, from 3.71 to 2.84 in the N 3 LO approx + NNLL approximation and from 2.95 to 
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2.85 in the N 3 LO approx ±NLO±NNLL one. The three-parameter fits strongly confine af’ (m|), 
within an error of 3.7% at 90% CL in both approximations. The inclusion of the r (Q 2 ) term has 
the beneficial effect of shifting ctf 5 (m|) closer to the world average, 0.1184 ±0.0007 [74]. In 
fact, our N 3 LO approx + NLO + NNLL value, 0.1199 ±0.0044 at 90% CL, which corresponds to 
0.1199 ±0.0026 at 68% CL, is in excellent agreement with the former. Note that similar (m|) 
valu has been otained recently [f76|| in an extension of the MLLA approach. 


7. Conclusions 


We have shown the ^-dependences of the SF F 2 at small-* values and of AJMs in the frame¬ 
work of perturbative QCD. We would like to stress that a good agreement wit the experimental data 
for the variables cannot be obtained without a proper consideration of the cotributions of both the 
“+” and ” components. 

The “+” components contain all large logarithms ln(l /x) as far as DIS SF F 2 and also for the 
average jet multiplicities. The large logarithms are resummed using famous BFKL approach [[jj] 
in the PDF case and another famous MLLA approach [50] in the FF case. 14 Nevertheless, the 
contributions of the ” components are very important to have a good agreement with experi¬ 
mental data: they come with the additional free parameters. Moreover, the ” components have 
other shapes to compare with the “±” ones. For example, in the AJM case the ” component is 
responsable for the difference in the (^-dependences of quark and gluon multiplicities. Indeed, the 
” component gives essential contribution to the quark AJM but not to the gluon one. 

In the case of DIS SF F 2 , our results are in very good agreement with precise HERA data at 
q 2 > 2 A 3 GeV 2 , where perturbative theory can be applicable. The application of the “frozen” and 
analytic coupling constants CCf r (Q 2 ) and a m ( Q 2 ) improves the agreement with the recent HERA 
data [|j] for small Q 2 values, Q 2 > 0.5 GeV 2 . 

Prior to our analysis in Ref. | [22| . [23| ], experimental data on the gluon and quark AJMs could 
not be simultaneously described in a satisfactory way mainly because the theoretical formalism 
failed to account for the difference in hadronic contents between gluon and quark jets, although 
the convergence of perturbation theory seemed to be well under control [25]. This problem was 
solved by including the ” components governed by 77 es (0, Q 2 -Qq) in Eqs. (5.10) and (5.12). 
This was done for the first time in Ref. j^2[]. The quark-singlet ” component comes with an 


arbitrary normalization and has a slow Q 2 dependence. Consequently, its numerical contribution 
may be approximately mimicked by a constant introduced to the average quark jet multiplicity as 
in Ref. @|. 

, 23] with fixed value of ocj 5 (m|), we then included 


ai 5 (m|) among the fit parameters, which yielded a further reduction of Xdo(- The fit results are 
listed in Table f|. 


Motivated by the goodness of our fits in | 


l4 Note, however, that in the case of DIS SF F 2 we use obly the first two orders of the perturbation theory and our 
“+” component resum by DGLAP equation [[]]. The resummation leads to the Bessel-like form of the “+” component. 
Including all orders of the perturbation theory should lead to a power-like form as it was predicted in the framework of 
BFKL approach [||j] (see discussion in [p4|]). 
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8. Appendix A 


For diagonalization of quark and gluon interaction it is neseccary to introduce the correspond¬ 
ing matrix U, which diagonalize exactly the LO AD 


U 


-1 ( Yqq{n) yfqin) 
W?(") Ygg\ n ) 


U = 


y-\n) 0 

0 yf\n) 


(U~ 1 U = 1), 


(Al) 


where 


U = 


a n a n — 1 

fin fin 


u~ l = 


1 -a„ 


_i is, 

. 1 A, 


(A2) 


and a n and fi„ have been defined in the main text, in Eq. (2.18). 

At higher orders the anomalous dimensions are transformed as follows 


IT 


yqq(n) 7gq(n ) 


7--O') 7-+W 
7?-(«) y++{n) 


if) 


where exact representations for y±±(n) and y± T (n) were given in the main text, in Eq. ( 2.22 ). 

8.1 Diagonalization of the renormalization group exponent 

Consider the renormalization group exponent (hereafter in the Appendix A a s = a s (Q 2 ) and 

a s = a s {Ql )) 


U = 


Si 


(A3) 


W{a s ,a s ) = T as exp 


C a s d a ' y[a') 
Ja s a' 2/3 («') 


(A4) 
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in the following form 


W(a s ,a s ) =V(a s ) exp 


r (0) W ln A 
2 A) ^s 


V-\a s ) = V(a s )W^(a s ,a s )V-\a s ), (A5) 


where the matrix V ( a s ) contains high order coefficients. 

To find the matrix V(a s ), it is better to find the derivation 


— W(a si a s ). 
da s 


The l.h.s. of (A5) leads to 


^-W(a s ,d s ) = ^^V(a s )W^(a s ,a s )V~ l (a s ). 


For the r.h.s. of (|A5|), we have 
d 


da 


■ W (a s ,a s ) = 


d Pn ^ t >t ^ 7 (0) (») 1 
— V(a s )-V(a s )—- - 

da s 2po a s 


Thus, the matrix V (a s ) obeys the following equation 

y(°) ( /2 


d ~, . 1 

1 ~V(a s ) + - 

da v 


2 A) 




JBT-P^W 11 V<*). 

2 p{a s ) 2p 0 a, / 


where the second term in the l.h.s. is the commutator of the matrices y ( °) ( n ) and V (a s ). 
Now we consider LO, NLO and NNLO approximataions separately. 

8.1.1 LO 

At LO, the matrix V (a s ) = I and the renormalization group exponent have the form 


where 


W(a s ,d s ) = W^{a s ,d s ) = 


d±(ri) = 



0 

_ \ d + {n) 
a s ' 


/ 


rl Q) W 

2/3 0 


8.1.2 NLO 

At NLO, the matrices V(a s ) and V ] (a s ) has the form 

V(a s ) =I + a s V^ l \ V~\a s )=I-a s V^\ 
and the Eq. ( |A9| ) can be replaced by one 

y(°) 


2t> (1) (n) + 


A) 


W(t 




A) 


A? 


(A6) 


(A7) 


# (0) (A8) 


(A9) 


(A 10) 


(All) 


(A12) 


(A13) 
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Applying the matrices U 1 and U to left and right sides of above equation, respectively, and 
using Eq. ( ]A3| ) for i = 1 and the representation 


U-'V^irDU = 


Vi'l {n) V%{n) 
V® + (n) 


for i = 1, we have the following matrix equation 

2v[ l l (n) vi!j* (n) • 2 ft+y- 0) w-y+ ) (« 

i-r^M 

'+-VV ft 


v ji)(„) 2ft+rl 0) (»)-r (0) M 


ft 

2 Vi%) 


J_(T ( -Hn) r«(n)' 

A) lr«(n) r«(«) 


where 


I±±0») = ri±(«) - 7± 0) (»)^1, r^(n) = 7±£(«) 


The Eq. (A 15) leads to the results 


Vfi (n) = vg(») = - 




2 A) + 7± ) («)-7t ) (») 


(A14) 


(A15) 


(A 16) 


(A17) 


8.1.3 NNLO 

At NNLO, the matrices V(a s ) and V 1 (a s ) has the form 

V(a s ) = I + a s V (1) + a 2 V {2 \ V~\a s ) = I- - a]V {2 \ y (2) = y( 2 ) - V {1) • V (1) (A18) 

and the Eq. (|A9|) can be replaced by one 


4t> (2) (n) + 


y(°) in) 

“AT 


1> (2) (n) 


A) L 


f- 2 \n) + y^ l \n) (v^^n) — (b\V^\n) +bj 


(A 19) 

Applying the matrices U 1 and U to left and right sides of above equation, respectively, and 
using Eqs. (A3) and ( [4 14 ) for i = 1 and i = 2, we have the following matrix equation 

4V®(n) }?[%). * ° + ^W-li 0, (") 

y (2) (n) .4M^Wz^M 4vi%) 


]_( r -- W r -+W 

A) \r;i(n) r$.(n) 


(A20) 


where 

p(2) 
1 ±± 


(«) = 7±±(«)+ E 7iV(«)^i 1) W-^i (Vii( n ) + 7± 0) ( n W±l ( ?1 )) - {b 2 -b\)y±\n), 

i=± v ' 


E±i(n) = 7±|M+ E 7i 1 i - ) W^ft 1) W-^ 1 ()©(«) + )i 0) («)V±=J(n)) 


(A21) 


The Eq. (A20) leads to the results 

rjiW t/(2) 


Vf2(n) = 


r£i(«) 


T/V Z 7 (yj\ — _ . 

4 A> ’ ±tU 4 j8 0 + 7± } in) — y£ 0) (n) 
42) 


(A22) 


Taking in brascets the relation between V and V^ 2 ) given in the last relation of (A18), 


we 


have 


Vi%) = V±±(n) - E Vg(n)v£>{n), V^[n) = Vjg(n) - E V^(n)V^(n), (A23) 


42), 




42) 


( 2 ), 


41)r„M/(l)/ 


i=± 


i=± 
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8.2 Q 2 evolution of parton distributions 

In the matrix form, the Q 2 evolution of parton distributions 

[f,(<2 2 ),f,(G 2 )] = lWl)MQl)\ -W(a s ,a s ) (A24) 

can be represented in the form 

[f,(e 2 ),f,(<2 2 )] = [WDMQo)] U • (U- x W{a Sl a s )U) U~ l . (A25) 

The first part in the r.h.s. is 

[fM)MQo)]U = [fM)^ + Wo)&MQo)(°^-l) + Wo)Pn] = [fs(<25),-fJ(G5)] , 

(A26) 


where (see also (2.17) in the main text) 


M<2o) = fq(Qo) a n+f g {Qo)Pn, ^(<2o) = Wo) (1 " «n) ~ Wo)Pn ■ (A27) 


8.2.1 LO 


At the LO, the renormalization group exponent U l W(a s ,a s )U has the diagonal form ( jA10| ) 
and, thus, we have 


[f M)MQl)\U ■ {U- l W(a s ,a s )U) = 


_ \ d-(ri) / — \ d+(n) 

C(<25)(-) -K(Ql) 


, (A28) 


Then, for the Q 2 evolution of parton distributions we have 


[W 2 )MQ 2 )\ = 


C(G5) 


2 \ / u s 


d-{ n ) 


-$(&) I 


d+(n)' 


1 TT 

1 


fa \ d -W fj; \ rf +(«) \-n fa \ d a fa 

\ M'S J \^S / Pn \^S / p/7 \ M'S J 




rf,(«) 


(=± 




2 \ / 


q 

diin) 


i=± 


where (see also (2.17) in the main text) 


r g {Q 2 ) =r q (Q 2 ) 1 —^ = i q (Ql)e n + i K (Q 2 ){\-a n ), 

Pn 

Pn 


because 


£„ = 


ex,,, (1 a n ) 
In 


(A29) 


(A30) 


(A31) 
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8.2.2 NLO 

At the NLO, the renormgroup exponent U~ l W(a s ,a s )U has the form 

U~ l W(a s ,a s )U = + - U^W^ia^a^U ■ (i-a s V^ 

Thus, it is convenient to consider firstly the part 

[f q (Q 2 o)MQo)]U-(l + a s V^ 1 


(A32) 


(A33) 


Following to ( |A26| ) we can rewrite it as 

[Wl)MQl)]U-(l + a s V^) = [f^(2o),-fJ(Go)] • (/ + a s F (1) ) = (g§) -^(Gg)] , 

(A34) 


where (see also (2.19) in the main text) 


if(».6o) = *f(«.Go) (n-OiVliW) ~^(n,Ql)a s V^l s {n) 


and 


V. 


T±,S 


— t/W y( 2 ) — 

- Fzp±, >^±,5 - k T± 


(A35) 


(A36) 


We introduce notations V^\_ s and Vwj: s in (A36), because the corresponding ones in the gluon 


case are different (see Eq.(A45) below). 
The Ov-part of (|A32) has the form 


U- l W (0 \a s ,a s )U (i—a. 


Thus, we have 




-a. 




-+ 




(A37) 


[Wl )MQl)\U-{U- l W(a s ,a s )U) 

' — \ d-{n) 

r(0 


77 \ d +( n ) 


t(Q§)(?) (i-«.v™)+a£(eS)(?) d 1 ? 


-{%«&) 


d+(n) 


i -d,v2) +d s r q {Ql) 


d -{n) 


V 


( 1 ) 


(A38) 


Then, to obtain the Q 2 evolution of parton distributions fs(<2 2 ) and f g(Q 2 ) we should product 
the r.h.s. of ( A38| ) on the matrix U 1 . By analogy with the calculations at LO, we have for quark 
density 


W 2 ) = rAQl) 


+^Ql) 


d-{n) 


1 -d s v[ V l ) + a#(Q; 


q \3<£.'o) 


d+{n) 


d+(n) 


1 -a s vi l i)+a s i-(Q z 0 ) 


V 


d -{n) 


(1) 


V. 


(1) 


-+ ■ 


(A39) 
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Taking together terms in the front of (a s /a s ) d± , we have 


r— \ rf_(«) 

W 2 )=t(Q 2 0 ){- 


or in more compact form 


1 — a. 


yd)-yd) 


+£(G§) 


'a \ <M«) 


1 - Cl, 


vi l l-v[ x i 


(A40) 


f q (n,Q 2 ) = £f'(«,<2 2 ), r q (n,Q 2 ) = f>,Go) -7 H' q (n,Q 2 ), 


i—zL 

H±(n,Q 2 ) = 1 -a. 


di(n) 


For gluon density we have 

f g (Q 2 ) = C(Go) 


1 ™ / — \ d_(rc) 

2 \ 1 ( a s 


Vl%)-V2l q (n) 

(l-a,V^)+5#(Go) 


(A41) 


l-On 

A? V ^.s / 


a /a \ d+ W 

2\ '-t'Yl I u s 


-t(Q 5)? (4) (1 -s.riS) - «A(e5)714 


-j2\ 


— \ d-(n ) 

n \ v ' 


q v u ' fi n \a s J V * ++ 7 

Taking together terms in the front of ( a s /a s ) d± , we have 

f g (<2 2 ) = C(Go) 


A? \®s, 


V 


( 1 ) 


-+ • 


(A42) 


i /— \ d-M 

1 — (Xn / a s ' 


-^(Go) 


Ai V d? 

n /a \ 

2\ / c \s 


1 — a* 


y(l) _ a » y(l) 

a„ -1 + 


fin \ ds 


1 -m 


or in more compact form 


vj$ - ——-vj 1 ] 

a„ 


77 \ 4 (n) 


(A43) 


f,(«,G 2 ) = £f>,G 2 ), f \{n,Q 2 ) = ?>,Go)HM //'(«, <2 


(=± 


^(«, G 2 ) = 1 - a s V£l ( n ) - (n), 


where 


r-+» = K-,M = 

and (see also ( 2.1 9| ) in the main text) 

ig(n,Ql) =i^(n,Qo)^-p^ L = f g"(«>Go) (l+a,vl‘2(n)) -£(n,Go) fl *^+-(») ^ 
= *7 («, Go) (i + («)) - (n, Ql)a s V+l (n) ——- 


r(l) / \ &vi 


rW 


- T_ 1 


(A44) 
(A45) 

(i)r„\ ^ ~ a " 


(Xn 


= *7 


(n,Qo) (l+fl,V^(n)) -£(GoM£>>), 


f^(«,Go) = ~^(n,Q 2 )^ =-f+(«,Go)| i (l + a,vj$9n0) + (n,Qo^vi? (n) — 


fin 


fin 


= £ {n, Go) (l + a,v|2(n)) - £ (n, Go)a*V-+ (n) 

= f)T(n, Go) (l + a.yj 1 ](n)) - (n, GoKvllj^ (n), 


(A46) 


or by analogy with ( |A35| ), in the general form, 

ff(«,Go) = ff(«,Go) (l+aA'EtW) - fj (n, Ql)a s V^l G (n) , (A47) 
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8.2.3 NNLO 


At the NNLO, we can perform an analysis, which is very similar to the one in the previous 
subsection for the NLO approfimation. The one difference is the terms ~ a 2 and ~ a 2 for the 


matrices V and V 1 (see Eq. ( A18 )). 

So, the final resuls have the form (a = q,g ) 


l,(n,Q 2 ) = £f>,<2 2 ), f>,<2 2 ) = f a (n,m ~ K(n,Q 


lb 


di{n) 


(A48) 


^(»>Qo) = tf(>hQl)(l+a s vi l l(n)+a 2 v¥±(n)^ -£(#i,Go) + «M±,«( n )) ■ 

8.3 <2 2 -dependence of Mellin moments 

The <2 2 -dependence of the singlet part M^{Q 2 ) of the Mellin moments can be obtained using 
the PDF Q 2 -dependence (see the previous subsection of the Appendix) and the relation between 
the parton densities and the (singlet part of) the Mellin moments given by Eqs. (IS) and ( 2.10 ). 

Sometimes, it is convenient to obtain directly the Q 2 -dependence of the Mellin moments 
Mf)Q 2 ). In the matrix form, it has the form ( A25|) 


(i) 


■A 2), 


>{ 2) 


C g (n,a s ) J 
(A49) 


M s n (Q 2 )=[f q (Q 2 ),f g (Q 2 )\ 

f Cq(n,a s ) 
yC g {n,a s ) 

where 

A-i ( QC 

vqo 

and 

C±{n,a s ) 

with 



C (ti. ci s ) 
-C+(n,a s ) 


(A50) 


(A51) 


B + («) = Bq *(») - tt B g ] (»), (' = 1 , 2 ) 

P/7 

B - (n) = B ^ (n) + 1—-^- B ( g ‘ ] (n ) = B^ (n) + (n ); 

P/7 P/7 


(A52) 


The basic idea is to split the (^-dependence to the initial conditions (Q)) ( |A27 ) and above 
LO to f^(«,<2o) ( |A35 ) and ( A48 ). The <2 2 -dependence combines the PDF one from the previous 
subsection and the one in (A51) and (A52). As it was above, we will consider LO, NLO and NNLO 
cases separately. 

Mote here that the Mellin moments M%(Qq) can be easy extracted from above equation (A49) 

MfM i) = [WeLMCo)] (?= HQl)MQl)](>■(>-' 


C g (n,a s ) 

= [r^(G5),-fJ(Ga)] = E f,(Ql)c,(n,c 


C g (n,a s ) 


(A53) 


i=± 


where Cj(n,a s ) are given by (A51). 
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8.3.1 LO 

Here 


M?(<2 2 ) = [f,(G 2 )4(<2 2 )] 


'-\d- in) 

£i ' 

. a s 

0 


0 

f- \ d + {n) 

£j ' 

. a. 



where 


8.3.2 NLO 


M^\Q 2 ) = \ l q (Ql)[ ^ 


d±(n) 


= (A54) 

i=± 


(A55) 


Here we have Eq. (A32) 

U~ l W(a s ,a s )U = (/ + a,y (1) ) •r 1 W (0 >(fl„fl s )E- (l-a s V (1) ) 


which leads to (J^-part (|A34j)) 


[f,(<26),f ? (<26)]^-(/ + ^ (1) ) = fc(G§),-i*(g§)], 


(A56) 


(A57) 


with ^(n, £>o) given by ( |A35[ ). 

The (> 2 -dcpcndcnt part consists from 


JJ- l W [0 \a s ,a s )U (l-a s V^ 


(A58) 


given by the r.h.s. of Eq. (|A37p and U l C(n,a s ) give by the Eqs. (|A5Q|) with the NLO coefficints 


(A51). 


So, we have 

U- l W {0 \a s ,a s )U (l-a s E (1) ) f/ _1 C(n,a. s ) 


({l )" W (l-^) 


V 

/ (aA d - {n) 


f- \d+{n) 


V 


(I) 


—a ' " 


77 \ d -(") 


( 1 ) 


V 

l-a,v{5 


l+a.B^ 


- 1 +a s B 


(i) N 


V U 


i t d s 


a \ d + ( «) , , 

s 1 + a s 


5 (i)_ya_) +y a_) 


(A59) 


For the Mellin moments M%(Q 2 ) we have product Q 2 and Q g parts: 


M n(Q 2 ) = I^ S,/) (<2 2 ), 


i=± 


where 


M^\Q 2 ) = f q {Ql)\ 


d±(n) 


1 +a s R i i ) 


(A60) 


(A61) 


39 

















Parton distributions 


Anatoly Kotikov 


and 


fi£ 1} = + 


Note that the eqs (A61)-(A62) can be rewriten as 


M^ ±] (Q 2 ) = 


'77 \ 1 + a s R 


? (i) 




1 +a s R { ± 


where 


M^\Ql) = ~^{Ql)(\+a s 


ff (G§) ( 


1 + Q.$ 




? (i) 


fi±+^1$ 


) +f5(2o) 


a s V^l 


8.3.3 NNLO 

Repeating the calculations in the previous case, we have 

(S- ! Wn 


MniQ 2 ) = Z^M^iO 2 ), 


i=± 


where 


,(i) 




rf±(;i) 


1 +a s R { 2 ] +a 2 R ( 2 ] 


'm = W) 

where f^(2o) an d are given by ( |A4S| ) and ( A62[ ), respectively, and 


>( 2 ) 


j(2) d( 1 ) f ..(I) 


R+ = fir - fl+ ( v±± - v±4) - vu + v. 


(i) 


>(2) , t/(2) 


r ±± T r± T 


A 2 ) 


t(2) 


with and V+j given by Eqs. ( |A23| ). 

The eqs (A66)-(A67|) can be rewriten as 


/ x d + (n) ( \ +d s R+ ) +d^R ( '+ ) ) 

A s >±)tn2\ - a s ’±)tn2\( a A V ^ ± ^ J ± ) 


M^\Q 2 ) = Mr ; (25)| 


X+a.R^+a 2 s R [ l ] 


where 


M, 


i s,±) (Qo )=?(G§)( 
ff(Go)( 


1 + a s R^± + a 2 R^± 


1 + Os 


B^+vQ 


+ a X 


fi^+fi^V^+yg 


+ fj (Go) 


(A62) 


(A63) 


(A64) 


(A65) 


(A66) 


(A67) 


(A68) 


(A69) 


di) 


a s V^+a 2 s vQ 
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